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Abstract 

We present details of earlier studies (Zheng et al, Phys. Rev. Lett. 78, 310 
(1997) and Das Sarma et al, ibid 79, 917 (1997)) and additional new results 
on double-layer quantum Hall systems at a total filling u = 2i^i, where a single 
layer at filling i^i forms a ferromagnetic, fully spin-polarized, gapped incom- 
pressible quantum Hall state. For the case vi = 1, a detailed Hartree-Fock 
analysis is carried out on a realistic, microscopic Hamiltonian. Apart from the 
state continuously connected to the ground state of two well separated layers, 
we find two double-layer quantum Hall phases: one with a finite interlayer 
antiferromagnetic spin ordering in the plane orthogonal to the applied field 
(the 'canted' state), and the other a spin singlet. The quantum transitions 
between the various quantum Hall states are continuous, and are signaled by 
the softening of collective intersubband spin density excitations. For the case 
of general vi, closely related results are obtained by a semi-phenomenological 
continuum quantum field theory description of the low-lying spin excitations 
using a non-linear sigma model. Because of its broken symmetry, the canted 
phase supports a linearly dispersing Goldstone mode and has a finite temper- 
ature Kosterlitz-Thouless transition. We present results on the form of the 
phase diagram, the magnitude of the canted order parameter, the collective 



excitation dispersions, the specific heat, the form of the dynamic hght scat- 
tering spectrum at finite temperature, and the Kosterhtz-Thouless critical 

temperature. Our findings are consistent with recent experimental results. 
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I. INTRODUCTION 



Interaction in a low- dimensional system does not merely result in strong renormalization 
of physical quantities, but can in many cases drive the system into completely new phases 
with peculiar properties. For a two-dimensional (2D) electron gas in a perpendicular mag- 
netic field, the interaction effects are especially important because of Landau level quantiza- 
tion. When electrons are entirely restricted to the lowest Landau level by a large magnetic 
field, electron-electron interaction completely dominates the properties of the system as the 
electron kinetic energy is quenched to an unimportant constant. One of the most interesting 
phenomena in this strongly-correlated system is the quantum Hall (QH) effect, which has 
attracted a great deal of experimental and theoretical interest during the last fifteen yearJ3. 
Recent advances in materials growth techniques have made it possible to fabricate high- 
quality double-layer two-dimensional electron systems with the electrons confined to two 
parallel planes separated by a distance comparable to that between electrons within a plane. 
With the introduction of this layer degree of freedom, many qualitatively new effects due 
entirely to interlayer correlations appeari^. Many new QH phases in double-layer systems 
become real possibilities because of the increased degree of freedom and the complicated 
interplay among interlayer tunneling energy, Zeeman energy, and electron-electron Coulomb 
interaction energy. 

In this paper, we present the details of our earlier theoretical investigationJBi of the 
possible QH phases in a double-layer system at a total Landau level filling factor u = 2z/i, 
where ui is a filling factor at which an isolated single layer system forms a fully spin polarized 
incompressible QH state {e.g. ui = 1, 1/3, etc.) We will discuss three distinct ground states, 
and the nature of the zero or finite temperature transitions/crossovers between them: 

• A fully polarized ferromagnetic (FPF) QH state in which the spins in each layer are 
aligned parallel to the magnetic field. This state is adiabatically connected to the 
ground state of well separated layers, each forming a polarized QH state at filling 
fraction z/i. We will denote this FPF state also as the FM (for "ferromagnetic") state. 
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• A spin singlet (SS) state, which can be visuahzed crudely as consisting of singlet pairs 
of electrons in opposite layers. Alternatively, at z^i = 1, we will discuss the Hartree- 
Fock picture of spin up and spin down electrons fully occupying single-particle states 
which are symmetric in the layer "pseudospin" index; hence the singlet state will also 
be referred to as SYM. In the limit of a vanishing tunneling matrix element between 
the layers, this state is simply the pseudospin polarized state of Refs for both 
spin up and spin down electrons separately. Throughout, we will consider the case 
of a no n- vanishing tunneling matrix element: in this case the pseudospin polarization 
is chosen by the phase of the tunneling amplitude, and not spontaneously. None of 
the phase transitions we consider here require a vanishing tunneling matrix element; 
on the contrary, changes in the value of the tunneling matrix element can drive the 
quantum transitions. 

• A canted state (C) in which the average spin moments in the layers have an an- 
tiferromagnetic correlation in the plane perpendicular to the magnetic field, and a 
ferromagnetic correlation parallel to the magnetic field. Both ferromagnetic and anti- 
ferromagnetic moments can vary continuously at zero temperature as parameters are 
varied. The planar antiferromagnetic ordering breaks spin rotation symmetry about 
the magnetic field axis: as a consequence there is a gapless, linearly dispersing. Gold- 
stone collective mode in its excitation spectrum and a Kosterlitz-Thouless transition 
at a finite temperature. The C phase is the canted antiferromagnetic phase (CAF) 
discussed in our earlier short publicationsi'i. 

We will use two distinct and complementary approaches to understand these phases. 
The first is a mean- field Hartree-Fock calculation: this applies only for integer values of z/i, 
but has the advantage of working with a precise microscopic Hamiltonian involving only 
parameters which are directly known experimentally. The second is a phenomenological, 
quantum field-theoretic formulation which applies for general z/i, and allows us to more 
precisely understand the consequences of thermal and quantum fluctuations. We will now 



discuss some of the results of these two approaches in turn. 

In the Hartree-Fock approximation0, we are able to show that the canted antiferromag- 
netic (C) phase is the energetically favored ground state ioi v = 2 at intermediate layer sep- 
arations for systems with small Zeeman energy, and that the phase transitions from the C to 
the FM or SYM phases are continuous. We evaluate at z/ = 2 the intersubband spin density 
wave (SDW) dispersions of all phases in the time- dependent Hartree-Fock approximation0 
and show that, as the precursor of the phase transitions, the collective intersubband SDW 
mode softens at the phase boundaries of the FM and SYM phases to the C phase. The 
SDW becomes the linearly dispersing Goldstone mode in the C phase, and the temperature 
of the Kosterlitz-Thouless transition is obtained by evaluating its effective spin-stiffness in 
the Hartree-Fock approximation. In addition we present results on the stability energetics of 
the various phases, the antiferromagnetic order parameter, the phase diagram, the collective 
intersubband SDW excitation dispersions, and the specific heat. 

The V = 2 Hartree-Fock results may also be qualitatively applicable to the case of z/ = 6 if 
the Landau level mixing is ignored (the Landau level mixing may not be negligible at u = 6, 
though.) On the other hand, the situation at = 4 is very different from the situation 
at u = 2, since the inter-Landau level excitation energies are comparable to the cyclotron 
energy; our results do not apply at z/ = 4. 

The microscopic Hartree-Fock analysis obviously does not apply to a situation where the 
average filling factor ui in each layer is fractional (e.g. ui = 1/3) with each isolated layer 
supporting a spin polarized Laughlin fractional QHE state; such a many-body state will not 
appear in any mean-field decoupling of the Hamiltonian. However, an essential property 
of the phases we are discussing is that they all have a gap towards charged excitations, 
and the transitions between them are driven by changes in the nature of the mean spin 
polarizations, and of the spin excitations. This suggests that it may be possible to develop 
a more general effective theory which focuses on the spin excitations alone. We will present 
such a theory in Section |ITT|: it turns out to be the 0(3) quantum non-linear sigma model in 
the presence of a magnetic field. For the case z/ = 2, we are able to use our earlier Hartree- 
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Fock computations to precisely obtain all the renormalized parameters which universally 
determine the low temperature properties of the non-linear sigma model; for other values 
of z/i, including the fractional cases, these parameters remain as phenomenological inputs. 
We will present the phase diagram of the sigma model, and describe the nature of the finite 
temperature crossovers above the various phases in some details. In particular, we will 
obtain explicit predictions for the temperature dependence of the line shape of the inelastic 
light scattering spectrum. 

We note that our findings from the two approaches are consistent with recent inelastic 
light scattering measurement!!, where a remarkable (and temperature dependent) softening 
of the long wavelength intersubband SDW mode in a z/ = 2 double-layer system is observed. 
We hope that our other explicit theoretical results may be tested in future experiments. The 



experimental situation will be discussed in Section |IV. 

This paper is organized as follows. The results of the Hartree-Fock theory are presented 
in Section In Section [11 A| , we study the ground state properties of the v = 2 double-layer 



system in a self-consistent mean-field approximation. The intersubband SDW excitations in 
the time-dependent Hartree-Fock approximation and associated mode softening are studied 
in Section [11 B| . The thermodynamic properties are discussed in Section [11 (J| , and some 



further discussion, along with an assessment of the validity of the calculation, appear in 



Sections [II D[ and [II E[ In a long and self-contained Section |T| we give our non-linear a 
model effective field theoretic description for a generic v = 2vi situation. Comparison of our 
theory with recent light scattering experiments is discussed in section |^ A short summary 
in section concludes this paper. We note that the readers who are interested only in 
microscopic Hartree-Fock theory could skip Section |T|, and the readers who are interested 
only in our long wavelength effective field theory could skip Section We have taken care 
in writing the two parts of our work, namely the microscopic Hartree-Fock calculation for 
V = 2 (Section |I|) and the non-linear a model description for v = 2vi (Section [1II[ ) as two 
separate self-contained pieces which can be read reasonably independent of each other if so 
desired. 
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II. HARTREE-FOCK THEORY 



We begin by writing down the explicit microscopic Hamiltonian of a double layer quan- 
tum Hall system. 

Within the lowest Landau level, the single particle eigenstates may be denoted by la/iO"), 
where a is the intra-Landau-level index in the lowest Landau level, /i = 0, 1 is the pseudospin 
index which labels the symmetric and antisymmetric subbands, and the spin index a = ±1 
labels t and | spins.lii The Hamiltonian of the double-layer system is 

H = Ho + H,, (2.1) 

where the non-interacting Hamiltonian is 

Ho = -Asas E( 1/2 - f^)Cl.C^,. - A, E ^Ci^.C^,^^ (2.2) 

where the pseudospin splitting Agas is the tunneling-induced symmetric-antisymmetric en- 
ergy separation, the spin splitting A^ is the Zeeman energy, and (C) is electron creation 
(annihilation) operator. The Coulomb interaction Hamiltonian Hi is 

cno2 ^ll^^2^l■i^^4, 0L10L2 q 

^C'^l+(2yAtl(TlC*Q2^2Cr2^a2+gHAt30-2C'ai^40-l, (2-3) 

where Vt is the area of the system, lo = {hc/eBy^"^ is the magnetic length. The non-zero 
Coulomb potential matrix elements are Vqooo = Vouo = ^looi = ^iiii = V+ and V^ioio = 
Vbioi = Viioo = K)oii = V_, with V±{q) = |K(g) ± Vbiq)], where Va{q) = ^ and Vbiq) = 
Va{q)e~'^'^ are the intralayer and interlayer Coulomb interaction potentials, respectively. (The 
finite well-thickness corrections can be taken into consideration by including appropriate 
form factors^.) 

The following subsections will examine various properties of if at z/ = 2 by mean-field 
and RPA-like treatments of the interactions in Hj. 
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A. Ground states 



In this subsection, we investigate the ground state properties of H, and obtain the three 
phases discussed in the Introduction. Performing Hartree-Fock pairing of (|2.3| ), one obtains 
the mean-field interaction Hamiltonian as 



where X,, 



-J— V VI/ 



5^ ^Ml/i20'10-2C*/'tl(Ti^M20-25 



(2.4) 

>, which depends on 



the electronic state being sought through the expectation value < Cj^,^^^C^^„-^ >. We self- 
consistently search for the symmetry broken states where, in addition to < Cj^^C^a >^ 0, 
the possibility that < C^|Ci_^| is also allowed. Because of the complete Landau 
level degeneracy, the Hartree-Fock Hamiltonian H^^ = Hq + H^^ in a uniform state is a 
4x4 matrix, representing the dimension of the subspace associated with the spin and layer 
degrees of freedom. It thus has four eigenenergies ei± and four eigenstates 0j± {i = 1,2), 
which are obtained as shown below. In the non-interacting base (|0 t))|l i))|0 T)); 
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(2.6) 



where 6'i and 62 are associated with the Hartree-Fock eigenstates (f)i± which need to 
be obtained self-consistently, ni± are electron occupation numbers {4>\±4>i±)-, and U± = 
^ X]p e~^^'°^^V±(p). The off-diagonal matrix elements Aj represent the possibility of the 
broken symmetry ((C*^|C*i-/xj,) 7^ 0) mentioned above. By diagonalizing the Hartree-Fock 
Hamiltonian H^^ of Eq. ( p.5| ), one obtains the eigenstates 

/ 



)1+,(PI-,(P2+,(P2- 



sin(0i/2) cos(ei/2) 
cos(^i/2) -sin(0i/2) 








and the eigenenergies 








E1 + E2 , 
E3 + Ei 
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sin(^2/2) cos(^2/2) 
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'{E,-E,y 



(2.7) 



(2.8) 



Eqns. ( p.5|) to (|2.^ ) form the complete self-consistent Hartree-Fock equations which need 
to be solved numerically. In fact, the only quantities to be determined in this self-consistent 
manner are the two parameters 61 and 62, which, in turn, uniquely define the eigenstates 
through Eq. (|2.7[). The eigenenergies always satisfy < ej+ = 1,2), so the ground 
state at = 2 is given by |) = Hj0|_|v), where |v) is the vacuum state. The ground state 
energy is given hj E = {Hq + ^H^^). 

There are several sets of 61 and 62 which make Eq. ( f^.T] ) the self-consistent solutions 
to the mean- field Hartree-Fock equations. One is ^1 = and 62 = 0, which corresponds 



9 



to the symmetric (SYM) state. Another is = and 62 = vr, which corresponds to the 
spin polarized ferromagnetic (FM) state. These two are the spin-ferromagnets (FM) or layer 
pseudospin-'ferromagnets'ilfl (SYM) whose existence is naturally expected in the presence 
of finite Zeeman and tunneling energies. More interesting is that, for Agas > ^z, there 
exists a solution at intermediate interlayer separations with < 6i < 71. As we shall see 
shortly, this new state possesses a canted antiferromagnetic ordering (the C phase), i.e. an 
interlayer inplane antiferromagnetic spin ordering with the inplane spin magnetic moment 
in each layer being equal in magnitude and opposite of each other. The energies of these 
different states are shown in Fig.|l|. It is clear from this figure that the energetically favored 
ground state is the SYM state at small interlayer separations, the C state at intermediate 
separations, and the FM state at large interlayer separations. The u = 2 double-layer QH 
system thus undergoes two quantum phase transitions as the layer separation is increased 
from (i = Oto(i— s>ooata fixed magnetic field. 

To show the antiferromagnetic spin correlations, we rearrange the eigenstates as 

<j,,± = {l/V2){\L)Sf^ + \R)Sl), (2.9) 

where s!^^^\ electron spin configurations in the left (right) layer in the eigenstate (j)i±, 
are ^^l = cos(^,/2)| t) - sin(^i/2)| j), Sf: = cos(ft,/2)| t) + sm{9i/2)\ i), and satisfy 
{Sl"_^ySlL = {Sf^ySfi = 0. We immediately obtain the canted antiferromagnetic spin order 

as 

{Sf) = -{S^} = ^(sin^i + sin^2), (2.10) 

where S^^^^ is the electron spin operator in the left (right) layer, and x denotes the spin 
alignment direction within the two dimensional plane. This canted interlayer antiferromag- 
netic spin ordering is shown schematically in Fig. ^ Note that the total spin magnetic 
moment still points in the direction of the magnetic field as required by symmetry. It is 
obvious that this antiferromagnetic order breaks the U{1) symmetry associated with the 
spin-rotational invariance of the system. Its consequences on the low temperature thermo- 
dynamic properties will be discussed later. The numerical result of this order parameter 
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l(iS^) — (iS^)l is shown in Fig.0. One can see that when Zeeman energy is increased, 
the range of the layer separations where the canted antiferromagnetic state exists shrinks in 
favor of the ferromagnetic state, as the Zeeman energy obviously favors the spin polarized 
state. It is clear that the phase transition is continuous. 

The phase diagram, shown in Fig. ^, can be constructed from this mean-field approx- 
imation. The states |0 t) and |1 ]) are occupied in the FM phase, |0 f) and |0 J,) are 
occupied in the SYM phase, and the C phase interpolates between them. The SYM phase 
exists for Asas > A^ and d < dd, the C phase exists for Agas > A^ and dd < d < dc2, and 
the FM phase exists for either A^ > Agas or > dc2- The FM phase is favored when A^ is 
increased, while the SYM phase is favored when Agas is increased. In the next subsection, 
the same phase diagram will be obtained by studying the softening of the intersubband SDW 
excitations in the time dependent Hartree-Fock approximation. 

In this subsection we have studied the ground state properties of u = 2 double-layer 
QH systems in a mean-field Hartree-Fock approximation and showed the existence of three 
stable QH phases. The most interesting observation is the existence of a canted antiferro- 
magnetic phase, with a broken spin rotation symmetry, in between the symmetric and the 
ferromagnetic phases. 

B. Intersubband SDW excitations and mode softening 

In this section, we study collective intersubband SDW spectrum of z/ = 2 double-layer 
QH systems in the time-dependent Hartree-Fock approximation.0 These excitations involve 
flipping both spin and pseudospin of the electron and are the lowest energy excitations at 
u = 2. The phase instability is studied by investigating the softening of the collective inter- 
subband SDW excitations. The results obtained in this section are in complete quantitative 
agreement with the results obtained from the ground state studies in the previous section, 
as, of course, they should be if the calculations are done correctly. 

In the absence of interaction, the two branches of the intersubband SDW excitations 
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which correspond to transitions |0 t) ^ |1 |) and |0 J,) ^ |1 |) have excitation energies 
|Asas ± where Ag^s and are interlayer tunnehng and Zeeman energies, respectively. 
Interaction renormahzes the excitation energies in two ways. One is a self-energy correction 
to the polarizability due to the loss of exchange energy when an electron is excited to a higher 
but empty level, which raises the excitation energies. The other is the vertex correction to 
the polarizability due to an excitonic attraction between the electron excited to the higher 
level and the hole it leaves behind, which lowers the excitation energies. In diagrammatic 
perturbation theories, the effect of the exchange energy on the excitation energies is ac- 
counted for by including the corresponding self-energy in electron Greens functions, and the 
effect of the excitonic attraction is represented by vertex corrections. The self-energy and 
the vertex correction must be consistent with each other obeying the Ward identity. The 
direct Hartree term does not influence the SDW excitations because Coulomb interaction 
is spin-rotationally invariant. Since the Coulomb interaction potentials are subband- index 
dependent, they may introduce mode-coupling between the two branches of the intersub- 
band SDW excitations. This mode-coupling pushes down the frequency of the low-lying 
excitation and hence helps mode softening. 

The intersubband SDW excitation spectra are obtained as the poles of the retarded 
intersubband spin-density response function 

X^<^\q,co) = -t e^-*([p5D(q,t),pL(-q,0)]), (2.11) 
Jo 

where the intersubband SDW operator is defined as 

P5D(r) = E0l_(r)0.+ (r). (2.12) 

i=l 

Psd{t') recovers to famihar forms psni^) = X]/^ C'^t(^)^i-Mi(^) ^Pi^ polarized state 

(^^1 = and 02 = tt), and Ps'd(i") = J2a C*oo-(i")C'i-cr(r) in the symmetric state (6'i = 62 = 0). 

X'''^*(g, u;) is evaluated in the time-dependent Hartree- Fock approximation ,0 which we 
adapt to double-layer systems and, for simplicity, we ignore all the higher Landau levels. As 
argued earlier, this should be a good approximation for our problem. In this approximation, 
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one includes the single-loop self-energy and the ladder vertex diagrams in the theory, which 
satisfies the Ward identities. This time-dependent Hartree-Fock approximation, therefore, 
corresponds to solving the vertex equation shown in Fig.^, where the electron propagators 
are the self-consistent Hartree-Fock Green's functions obtained from the mean-field approx- 
imation discussed in the previous section. Due to the fact that the Coulomb interaction is 
frequency independent and that the Landau levels are completely degenerate, the integral 
vertex equation can be transformed into an algebraic matrix equation.0 The matrices can 
be further block diagonalized into 4x4 matrices, from which the poles of the spin-density 
response function can be (almost) analytically calculated. 

Combining Eqns. (|2.11| ) and ( p.l2| ), one obtains the spin density response function in 
the Matsubara frequencies0 

X{q,too) = e-^^'-/^ Ve-^«-°'°A+(^c.)r,+ (g,2u;,a), (2.13) 



where 



P n 

iuj + Ei-x - EiX 

= — i for T = 0, (2.14) 

AlUJ + Ei^ — Ei+ 

where j3 = l/ksT, Qix is the Green's function corresponding to the self-consistent Hartree- 
Fock eigenstate (pix and eigenenergy Eix given in Eqns. ( p.7|) and ( |2.8| ), respectively. The 
ladder diagram vertex function is 

riA(g,2u;,a)= e''?-"'° - ^ e-\vl+{^-^'?\ill'^^P.ivil (2.15) 
xDi>x'Ti>x'{q,iuJ,a'){iX;i' - X'\V{p^,a- a')\i - X]i'X'), 



where the interaction matrix element is 



(^iAi;z2A2|^(g)K3A3;a4)=^[l + (-l)'^+^^+-+-] ((^.^aJ^^.-aJ ((^.- aJ^^.- A3 



2' 

x[VMK^^+V-{q){l-S^.^s)], (2-16) 
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where Sji^ is the electron spin states given in Eq. ( p.9|) . 

To solve the vertex equation, we perform the following Fourier transformationsS 

Ta{k) = J2Ta{a)e-''^'o^ (2.17) 

a 

and 

V^.V'A'(q) = ^ E e-^''°/'e^P^'''^V^,,,,v(p), (2.18) 
" p 



where p A q = p^Qy - PyQx and ViX;i'x' = {iX; i' - X'\V{q)\i - A; i'A'), as given by Eq. ( p. 161) . 
After an analytical continuation, one obtains 

X^^*(g,^)=e-^''^/2 E T.+ (g,^), (2.19) 

i=l,2 

where 

T = (D-^ + y)-^Ar, (2.20) 

and T are 4x1 matrices, with Nix = Q/ 27:11, number of magnetic flux passing 
through the system, and Ta = Dix{uj)Tix{q,uj). D and V are 4x4 matrices, with Dix^i'x' = 
6ii>6xx'Dix{uj), and Vixyy deflned in Eq.. ( |2.18| ). 

The intersubband SDW dispersion uj{q), which occurs as the pole of the retarded spin 
density response function x^'^^^ is the solution to det\D^^{uj) + V{q, uj)\ =0. After a lengthy 
but straightforward algebraic manipulation, the two intersubband SDW dispersions uj±{q) 
are obtained as 

Lul = A^ + B^~ cos(0i + 62) 



± 



2 



{1 - cosiOi + 62)) A +4B^{A + C){A-C) , (2.21) 



where A = ^{a + b), B = \{a - b), and C = (1 + cos(6'i + 6'2)), with 

a = y^(Asas + A, + U+ cos Oi - f/_ cos 62)^ + {U+ sin Oi + f/_ sin 62)^ - V+, 

b = ^/(Asas - A, + f/+ cos 02 - f/- cos ^1)2 + (f/+ sin 62 + U. sin O^y - V+. (2.22) 
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The intersubband SDW dispersions in both the canted antiferromagnetic QH phase (C) 
and the normal QH phases (FM or SYM) can be obtained from the above expression by 
incorporating appropriate values of di and 62. In the following, we show uj±{q) only at 
zero temperature for the sake of simplicity, although the formalism applies equally at finite 
temperatures. 

In Fig. H, we show the dispersion of the intersubband SDW above the FM ground state. 
As mentioned earlier, these two intersubband SDW modes u!±{q) correspond respectively to 
transitions |0 t) |1 i) and |1 t) — |0 |). The frequencies u!± increase as functions of q, 
approaching asymptotic values uj±{q — > 00) = 00^ + \vx\, where are the non- interacting 
excitation energies and is the exchange energy of the electron in the ground state. Mode 
coupling, which pushes down uj^{q) and hence helps mode softening, is most visible at 
g — ^ 0. At zero layer separation, mode-coupling disappears, and we recover previously 
known results.0lll In Fig. |^, we show the intersubband SDW dispersion above the SYM 
state. The results are qualitatively similar to those in Fig. ^, except that there is no mode 
coupling in the symmetric state because Coulomb interaction is spin independent. The 
important thing to be noticed is that the long wavelength collective excitations are gapped 
in both the symmetric phase and the spin polarized phase. However, the mode softening 
does occur at the phase boundaries, as we show below. 

To illustrate the phase instability, we show, in Fig.§, the lower-energy branch of the 
intersubband SDWs at g = as a function of interlayer tunneling. We see that uj-{q = 0) 
indeed softens when approaching the phase boundaries from both the symmetric phase and 
the spin polarized phase, and remains zero inside the canted antiferromagnetic phase. The 
canted antiferromagnetic order parameter, calculated in the previous section, is also shown 
in Fig.^ for comparison purpose. We notice that the phase boundaries determined from 
these two independent approaches agree completely, as shown in the figure. The softening 
of the collective mode and the appearance of the antiferromagnetic order parameter implies 
that we have discovered a quantum phase transition in double-layer QH systems. 

In Fig.p, the collective intersubband SDW dispersions in the canted antiferromagnetic 
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QH state are shown. The first thing to be noticed is that the lower energy branch uj-{q) is a 
gapless mode. The existence of such a gapless Goldstone mode is due directly to the canted 
antiferromagnetic spin ordering which spontaneously breaks the spin-rotational symmetry 
of the Hamiltonian. This Goldstone mode is found to be linear in the long wavelength limit, 
consistent with the fact that it describes antiferromagnetic fiuctuations. The existence of the 
gapless excitation in the canted antiferromagnetic phase implies that some thermodynamic 
quantities, such as specific heat, have power-law temperature dependence in the canted 
antiferromagnetic phase in contrast to their exponential temperature dependence in the 
normal (symmetric or ferromagnetic) phases. 

Simple expressions governing the phase boundaries can be derived from the mode soften- 
ing. The boundary between the symmetric (SYM) phase and the canted antiferromagnetic 
(C) phase is found to satisfy the following equation 

(A3as-f/-)' = f/' +A^, (2.23) 

where f/_ = V^{q = 0) = ^ 6^*'^'°^^V"_(p). It should be noted that, for any given Agas, 
the critical layer separation at this boundary is considerably smaller than the critical layer 
separation where the charge density excitation in the v = 1 state becomes soft.lii The reason 
for this is the absence of Hartree contribution to the SDW excitations. The boundary 
between the spin polarized (FM) phase and the canted antiferromagnetic (C) phase is found 
to satisfy 

(A, + [/_)2 = [/!+AL. (2.24) 



The simplicity of Eqns. ( |2.23| ) and (|2.24| ) makes the phase diagram easy to construct. It 



is worthwhile to note that the phase boundaries are determined by only three energy scales 
A^, Asas, and f/_ in spite of the fact that the Hamiltonian is determined by four independent 
energy scales A^, Agas, and V±(g), of which the inter- and intra- layer interactions V±{q) are in 
fact continuous functions of wavelength q. This unexpected dependence of the phase diagram 



(Fig. |iy) on just three energy scales which are entirely determined by the magnetic field, the 
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sample parameters (i.e. inter-layer separation, well width, etc.), and the tunneling strength, 
is a specific result of the Hartree-Fock approximation. The zero temperature phase diagram 
can thus be expressed as a function of two independent dimensionless variables A^/Asas and 



f/_/Asas, as shown in Fig. |T0|. This phase diagram applies to all double-layer quantum Hall 
systems at z/ = 2 which may have any values of Zeeman energy, tunneling energy, layer 
separation, layer-thickness, etc. We believe, however, that this remarkable scaling in the 
phase diagram (which enables us to reduce an infinite number of Asas versus d diagrams 
for various values of A^, of which examples are shown in Fig. to just one phase diagram 
shown in Fig. [T^) remains approximately valid, although the relative size of various phases 
in the universal phase diagram of Fig. ^ may very well be quantitatively not particularly 
accurate. We also mention here that this phase diagram is topologically identical to that of 
a (2 + 1) -dimensional quantum 0(3) nonlinear cr-model in a magnetic field,! as discussed in 
section |I| of this paper. 

In this subsection, we have studied the collective intersubband SDW excitations for 
u = 2 double-layer QH systems in the time dependent Hartree-Fock approximation. We 
have presented the dispersions of the collective SDW excitations in both the normal QH 
phases (FM and SYM) and in the canted antiferromagnetic QH phase, and investigated the 
mode softening which signals the phase instabilities. We have rederived the same phase 
diagram as that obtained in the previous section, and obtained analytic equations for the 
two phase boundaries separating the new canted antiferromagnetic phase from the normal 
FM and SYM phases. 

C. Kosterlitz-Thouless transition 

In this subsection, we discuss some thermodynamic properties of z/ = 2 double-layer 
systems which arise from the spontaneous symmetry-breaking associated with the breaking 
of U{1) planar spin rotational symmetry in the canted antiferromagnetic quantum Hall 
phase. There should be a finite temperature Kosterlitz-Thouless transition in the canted 
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antiferromagnetic phase, since the spin-rotational U{1) symmetry is broken. Below the 
critical temperature, the system supports a linear Goldstone mode, which gives rise to a 
power-law temperature dependence for the specific heat. Above the critical temperature 
the U{1) symmetry is restored and the system is paramagnetic. These properties are, in 
principle, experimentally observable and provide direct ways to test our theory. 

We can estimate the Kosterlitz-Thouless transition temperature for our problem in the 
following manner. In the canted antiferromagnetic phase, the low temperature thermo- 
dynamics is governed by long wavelength phase fiuctuations of the order parameter. Let 
Eff) = {(j)\H\(f)) — where |) is the ground state of the canted antiferromagnetic phase, 

and 10) = exp{iJ2jSj4>j)\), with as the spin operator of the j-th electron and z is the 
(magnetic field) direction normal to the two dimensional plane. In the long wavelength limit, 
one obtains 



Ic/Vl v0(r)|', (2.25) 



with 



OO on. 



, , / sin 6'i + sin 6^2 \ ^ , . f sin 6i — sin 6*2 \ ^ 
Vaiq)[^ 2 ) +^^(^)( 2 ) 



dq, (2.26) 



where lo is the magnetic length and Va (fe) is intralayer (interlayer) Coulomb potential. 
For future convenience, we have written the stiffness as an explicit function of A^, which 
arises from the dependence of the angles 61^2 on the Zeeman splitting. The effective pla- 
nar XY model defined by Eq. (|2.25| undergoes a Kosterlitz-Thouless phase transitionB at 
approximately /c^^c = {'^/'2)Ps{^z)- Finite temperature spin- wave and vortex-antivortex 
polarizations reduce the transition temperature to approximatelyi'i0 

keT, ^ 0.90p,(A,). (2.27) 

These finite temperature renormalizations can be much larger in the vicinity of the C-N 
and C-SYM phase boundaries: the expression (^.27]) can then no longer be used, and we 



will discuss modifications near these boundaries later in Section III 



Since we know ^^(A^) exactly within the microscopic Hartree-Fock approximation, the 
Kosterhtz-Thouless transition temperature can be easily determined for our problem. In 
Fig. |lT], we show the calculated Kosterlitz-Thouless critical temperature in u = 2 double- 
layer quantum Hall systems within the mean-field Hartree-Fock approximation (i.e. from 
Eq. ( |2.26D for ps{Az)). The phase transition exists only in the canted antiferromagnetic 



quantum Hall phase. The critical temperature vanishes at the phase boundaries as the 
symmetry-breaking order parameter drops continuously to zero as the phase boundaries are 
approached from within the canted antiferromagnetic phase. We notice that the calculated 
Kosterlitz-Thouless temperature (~ IK) is well within the experimentally accessible regime 
for typical AlGaAs/GaAs-based double-layer systems. The effective spin-stiffness ps{Az) 
given in Eq. ( p.26|) is obtained in the mean-field Hartree-Fock approximation, i.e. using the 



results from sections |II A| and [II B|, where quantum fiuctuation effects are not included. The 



results in Fig. |TT| should thus be regarded as the upper bound for the Kosterlitz-Thouless 
critical temperature. We emphasize that the Kosterlitz-Thouless transition discussed here is 
present even in the presence of interlayer tunneling (in fact, the presence of finite interlayer 
tunneling is essential to stabilize the canted antiferromagnetic phase, as described in the 
last two sections), unlike the case associated with the pseudospin transitionifi at u = 1/m 
(m odd integers) where interlayer tunneling suppresses Kosterlitz-Thouless transition. 

Below the Kosterlitz-Thouless transition temperature, the specific heat in the antiferro- 
magnetic phase has qualitatively different temperature dependence from those of the normal 
quantum Hall phases. This is of practical significance since it is possible to experimentally 
measure the specific heat of quantum Hall systems.^ At low temperatures, the main con- 
tribution to the specific heat comes from long wavelength low energy intersubband SDWs. 
With their dispersions calculated in each of the quantum Hall phases, the heat capacity is 
easily obtained: C = {d/dT)J2k{n-{k))u!^{k), where uj^{k) is the energy of the low-lying 
intersubband SDW excitation and {n^{k)) is its Bose occupation factor. The results are 
shown in Fig. |T2[ It is clear that the specific heat has an activated behavior in the normal 
quantum Hall phases because of the existence of excitation gap in its spin wave spectra, 
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and a quadratic power-law temperature dependence in the canted antiferromagnetic phase 
because of the existence of the hnear Goldstone mode in the symmetry broken phase. The 
spin stiffness goes to zero discontinuously at Tc, and for T > Tc we have the usual disordered 
X-Y phase of the Kosterlitz-Thouless transition. 



D. Multicritical points 

Our analysis so far has obtained solutions for the FM, SYM and C phases obtained by 
varying the parameters A^, Ag^s, d in the Hamiltonian (see Figs. | and|l^) which modify the 
relative strengths of the Zeeman energy, the tunneling energy, and the Coulomb interaction 
energy, respectively. Generically, these phases separated by phase boundaries representing 
second-order quantum transitions. However, there are also special quantum multicritical 



points in Figs ^ and |iy whose physical significance we will now discuss. 

The first quantum multicritical point is apparent in Fig. ^ where the FM, C and SYM 
phases come together at a single point on the abscissa corresponding to vanishing interlayer 
separation (d = 0). It is easily noted from Fig. |^ that this quantum multicritical point is in 
fact defined by 

Asas = A,; d = 0, (2.28) 
which is equivalent to the conditions 

Asas = A,; V^{q) = 0, (2.29) 



using the definition of V±{q) given immediately following Eq. (|2.3| ) in section [ll7\| . The 



simple physical reason for the vanishing of V- (q) along the d = line is that the intra- and 
interlayer Coulomb interactions are identical in the limit of vanishing interlayer separation 
d. Note also that the vanishing of V- ( and consequently of f/_) pushes the quantum mul- 
ticritical point to an infinite value of the abscissa ( Agas/f^- oo) in the scaled universal 
diagram given in Fig. |TU] — in Fig. |TU| the two phase boundaries separating the three phases 
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approach each other asymptotically as Asas/f^- — ^ oo and A^/Asas — > 1 at the multicritical 
point. Note that the condition A^ = Agas for the quantum multicritical point is a particu- 
larly interesting criterion because, in the absence of our predicted canted antiferromagnetic 
phase (i.e. if the u = 2 double-layer QH systems allowed only the ferromagnetic and the 
symmetric phases, as was assumed in the literature before our work), the condition of the 
equality of the Zeeman splitting and the symmetric-antisymmetric gap (i.e. A^ = Asas) is 
precisely the single particle level crossing criterion where, at = 2, one would make a tran- 
sition from the ferromagnetic phases where the two up-spin symmetric and antisymmetric 
levels are occupied and the down-spin levels are empty for A^ > Agas to the symmetric (spin 
singlet) phase where the spin-up and spin-down symmetric subbands are occupied (and the 
antisymmetric levels are empty ) for Agas > A^. What our theory definitely predicts is that 
such a simplistic one particle level crossing picture (which appears to be obvious intuitively) 
does not occur in a double-layer QH system at z/ = 2 — instead Coulomb interaction breaks 
the SU{2) spin rotational symmetry and drives the system into an intervening antiferro- 
magnetic phase where spin and pseudospin levels are intrinsically mixed. The fact that the 
intuitively expected level crossing phenomenon (at A^ = Agas) has never been observedS 
m. a. u = 2 double-layer QH system in spite of systematic effortli is, in our opinion, rather 
strong indirect evidence in support of our phase diagram. 

The second multicritical point becomes apparent only in the universal phase diagram 
shown in Fig. [1^ (and can be inferred implicitly from the trend that can be seen in the 
phase diagrams shown in Fig. ^). Its existence is a consequence of the intriguing finding 
that our antiferromagnetic state, in fact, persists all the way to A^ = (as can be clearly seen 



in Fig. |T0| where a finite region of the antiferromagnetic state exists along the A^ = line) 
where the spin-polarized ferromagnetic phase no longer exists, and the antiferromagnetic 
phase is separated from the spin singlet phase by a multicritical point (M) defined by the 
condition 

Agag = 2U- with A^ = 0. (2.30) 
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Thus the critical hne defining the phase boundary between the antiferromagnetic and the 
symmetric phases for 7^ ends at a critical point (M) for = 0. It is evident that in 
the absence of any Zeeman energy (A^ = 0) the spin magnetic moment in each layer lies 
completely in the 2D plane of the electron gas where they must be equal and opposite in 
the two layers. Therefore, the A^ = antiferromagnetic phase of Fig. |TD| is not a canted 
phase, but is a purely Neel phase (N); indeed the Hamiltonian has full SU(2) spin rotation 
symmetry for A^ = 0, and spin moments in the N phase can point in any two anti-parallel 
directions. The N, C, and SYM phases meet at the multicritical point M. This multicritical 
point will take on special significance in our effective field theoretical formulation in the next 
section. 

Let us also note that the existence of this purely Neel QH antiferromagnet at z/ = 2 
double-layer system may not be just a theoretical curiosity because it is possible to obtain 
vanishing Zeeman splitting in a GaAs double-layer system in a finite magnetic field situation 
by applying external pressure which under suitable conditions could lead to the vanishing 
of the effective gyromagnetic ratio (the (yf-factor) due to band structure effects. 

E. Comparison with earlier work 

Before concluding this section, and going on to the effective field theoretic description of 
the double-layer QH system, we will discuss the relationship of our results to some earlier 
work on double layer systems. We will also use this opportunity to comment on the validity 
of the Hartree Fock approximation in our and earlier work. 

Most earlier studies,i^'@ however, have focused on u = 1 (with some work! on u = 
1/2). Although the u = 1 and the u = 2 QH systems exhibit some similarities such as 
the softening of their low energy collective excitations under certain conditions, there are 
important distinctions between them. At u = 1, the spin degree of freedom is normally frozen 
out by the external magnetic field. The relevant low energy excitations in the u = 1 QH 
state are therefore intersubband charge-density-wave excitations, and the properties of the 
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system are determined by the interplay between interlayer tunneling energy and Coulomb 
interaction energy. In this sense, the u = 1 system is in fact a single-layer system with a layer 
pseudospin-dependent interactionSBEl At z/ = 2, both the spin degree of freedom and the 
layer degree of freedom are relevant, and the low energy excitations are intersubband SDW 
excitations. Consequently, the properties of the system are determined by the interplay 
among tunneling energy, Zeeman energy, and Coulomb interaction energy. Because of the 
increased degree of freedom, the system has more ways to optimize the total energy, and new 
states which are not possible at u = 1 become possible at u = 2. The symmetric QH state 
is energetically favored at small layer separations because it optimizes the tunneling energy; 
The spin polarized QH state is favored at large layer separations because it optimizes the 
Coulomb interaction energy; The canted antiferromagnetic state is energetically favored at 
intermediate layer separations. The reason for this is that the canted antiferromagnetic state 
tends to simultaneously optimize both the tunneling energy and the Coulomb interaction 
energy, which prevails at intermediate layer separations where the tunneling energy and the 
Coulomb interaction energy are equally important. Both the canted antiferromagnetic state 
and the symmetric state exist only for systems with small enough Zeeman energy, as the 
Zeeman energy clearly favors the spin polarized state. 

Another important distinction between the u = 1 systems and the u = 2 systems is that 
although at u = 1 the mode softening destroys the QH effect^, and there is no reliable 
description of the electronic state in the non-QH phase because beyond the critical layer 
separation the system becomes effectively a pair of isolated layers with compressible half- 
filled Landau level states; in contrast, at z/ = 2, the QH effect prevails at all phases because 
there is always a charge gap in both layers (even as d ^ oo), and we have good understanding 
of the ground state and the low energy excitations in each phase due to the existence of 
incompressible filled Landau levels. Nevertheless, the mode softening and the associated 
phase transitions at u = 2 are likely to be observable through inelastic light scattering 
experiment^Uii and thermodynamic measurements. 

Our work has studied u = 2 double-layer systems by numerically solving the self- 
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consistent mean-field equationscj, and obtained collective excitation dispersions using many- 
body diagrammatic techniquesJlil Both approaches are, however, based on the Hartree-Fock 
approximation. In single-layer integer QH systems, calculationJi^ in the Hartree-Fock ap- 
proximation agree well with experiments. In double-layer systems, the Hartree-Fock ap- 
proximation is less accurate because Coulomb interaction potential is layer-index dependent. 
Nevertheless, we expect that the Hartree-Fock approximation remains a reasonably good de- 
scription for a double-layer system at z/ = 2, since the Hartree-Fock ground state, which is 
non-degenerate and separated in energy from higher levels, is a good approximation for the 
real many-body ground state at z/ = 2 due to the existence of incompressible filled Landau 
level states with charge excitation gaps at any layer separations. We want to especially 
emphasize the difference in the validity of the Hartree-Fock approximation between u = 1 
and u = 2. The approximation is valid at u = 1 only at small layer separations and fails 
completely beyond a critical layer separation where the system becomes effectively a pair 
of isolated layers with compressible half-filled Landau level state in each layer. At = 2, 
incompressible states with filled Landau levels exist at any layer separations. In particular, 
there is still one filled Landau level in each layer at c/ oo. This fact, namely the exis- 
tence of an incompressible energy gap at all layer separations, ensures that the Hartree-Fock 
approximation, upon which our calculations are based, is a reasonable formalism at z/ = 2 
regardless of the value of the layer separation. 



III. CONTINUUM FIELD THEORY AND QUANTUM CRITICAL PHENOMENA 

The Hartree-Fock analysis used in the previous sections has the advantage of working 
with a realistic microscopic Hamiltonian and of making definite quantitative predictions for 
experimental observables in realistic samples. In this section, we will present an alternative 
analysis based upon a continuum effective quantum field theory for the low-lying spin exci- 
tations of a double layer quantum Hall system. We will find that the global phase diagrams 
obtained in the two approaches are very similar, and are, in fact, topologically identical, and 
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that detailed additional predictions for the temperature dependence of various observables 
can be made by a combination of the two methods. In particular, some advantages offered 
by the continuum approach are: 

• It will become clear from the analysis below that there are two basic ingredients 
necessary to obtain the phases of Fig [1^: two well separated layers form fully polarized 
ferromagnets with a gap towards charged excitations (i.e. an incompressible QH effect 
gap), and the primary coupling between the layers is an antiferromagnetic exchange 
(i.e. a superexchange) interaction. As such, we expect a similar phase diagram to 
apply not only at filling u = 2, but also at any u = 2z/i, where z/i is any filling fraction 
where a single layer has a charge gap, and is fully polarized. In particular, this criterion 
is satisfied at = 2/m, m an odd integer, where each layer forms a polarized Laughlin 
fractional quantum Hall state. The Hartree-Fock analysis clearly cannot be applied 
for m > 1 , as the single layer charge gap appears only after inclusion of the non-trivial 
correlations implicit in the Laughlin state. 

• The Hartree-Fock theory significantly overestimates the energy of the spin-unpolarized 
symmetric (SYM) or the spin singlet (SS) state, as we will refer to it in this section. 
Spin up and spin down electrons are simply placed into the same orbitals which are 
symmetric in the layer index. This is costly in Coulomb energy as there are no cor- 
relations in the layer positions of the up and down spin electrons. It is clearly more 
advantageous to form a spin singlet states between electrons which are localized in 
opposite layers. The non-linear sigma model continuum field theory to be discussed 
below does this in a natural way. From now on in this section we refer to this symmetric 
or the spin singlet phase as the SS to emphasize its correlated singlet nature. 

• A number of quantum-critical points have been uncovered in the Hartree-Fock analy- 
sis. There is the = quantum critical point between the spin singlet (SS or SYM) 
and the Neel (N) phases, and a critical line between the SS and the canted antiferro- 
magnetic (C) phases. Our continuum approach will obtain the critical theory for these 

25 



transitions, and we will find that they have dynamic critical exponents^ z = 1 and 
z = 2 respectively. There is also a second critical line between the C and the fully spin 
polarized ferromagnetic (FM) phases: this transition has z = 2 and will be discussed 
only in passing, as the critical theory is rather similar to one of the models discussed 
in detail in Ref 

• The continuum theory offers not only provides us the zero temperature quantum phase 
diagram but also a streamlined approach to the study of properties at non-zero tem- 
perature, especially in the vicinity of the quantum critical points where effects of 
fluctuations cannot be neglected. The price one pays is that in general the param- 
eters defining the effective field theory are quantitatively unknown and can only be 
calculated from a microscopic theory such as the Hartree-Fock theory of the previous 
sections. 

We motivate our formulation of the continuum theory by consideration of the physics 
of two well-separated identical layers at u = 2/m. More specifically, the layer separation, 
d, is much larger than the magnetic length, io- Then the two layers (labeled 1,2) are 
approximately decoupled, and each separately has filling fraction ui = V2 = ^/fn. Their 
ground states will be the familiar Laughlin states for m > 1, or a fully filled lowest Landau 
level at m = 1, both of which are incompressible states with large energy gap to all charged 
excitations. These states are also fully spin polarized; the spin polarization is induced not 
just by the Zeeman coupling to the external magnetic field, but also by the significantly larger 
intralayer ferromagnetic exchangeiliS. The low-lying excitations in each layer are spin waves 
which have a small excitation gap given precisely by the Zeeman energy = gfisB. For 
small g, a complete descriptionllil of the low energy excitations of each layer can be given 
in terms of an action for unit vector fields ni^2 {^12 — ^) representing the orientation of the 
ferromagnetic orders. Spin waves are small fluctuations of ni,2 about an ordered state, while 
charged quasiparticles are SkyrmionBSHi textures of ni 2. The effective action describing 
the two layers isi'Silii (in units with h = = 1) 
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Lp\n\ = iMoA{n) ■ |^ + y (V.n)' - MoA.n, (3.1) 



Here 



'''' = 4^ 

is the magnetization density per layer, with Iq the magnetic length. The spin stiffness 
of each well separated layer is represented by p°; for m = 1, we have the exact result0 
Ps = e^/ (16V2vre£o), while for m > 1 numerical estimates of p° are given in Refil. The term 
involving A accounts for the Berry phase accumulated under time evolution of the spins; 
here A is any functional of n which satisfies 

dA,{n) 

6,-.^^ = n. (3.3) 

This Berry phase term also has a "dual" interpretation in the picture in which Cp is viewed 
as an action for Skyrmionsffl: it represents the coupling of the Skyrmion current to a 
"magnetic field" of strength AirMo. 

Now imagine reducing the value of d to couple the two layers. As there is a charge gap in 
each layer, we can neglect all charge transfer processes, and focus solely on spin exchange. 
Because of the strong repulsive interactions within each layer, we expect by an extension 
of the familiar arguments made in the context of the Hubbard model that there will be 
an antiferromagnetic superexchange coupling between the layers. This can also be inferred 
easily by considering the leading effect of interlayer tunneling and Pauli principle, which 
immediately provides a superexchange coupling between the layers. The complete double 
layer action is therefore 

Si = J d^x J dr {Cplni] + CF['n2] + Jni ■ 712) (3.4) 

The value of the interlayer exchange, J, is not known precisely; we expect that it is of order 
J ~ MqA1^/U where Asas is the tunneling matrix element (see Eq. ( p. 2]) , for example) 
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between the layers, and U ~ e'^/eio is the Coulomb repulsion energy. In addition to the 
imprecisely known J, the present approach also requires knowledge of the nature of the short 
distance cutoff at lengths of order io beyond which present continuum approach cannot be 
applied. We will show later that our ignorance of these quantities can be reduced entirely 
to uncertainties in the value of a certain renormalized energy scale. This energy scale can 
be either measured directly in an experiment, or computed by any microscopic theory such 
as the Hartree-Fock approach (appropriate at z/ = 2) described in sections [II AHIIC^ . Apart, 
from this single energy scale, however, all of the predictions of the present effective field 
theoretical approach will be quantitative and precise. 

Some potentially important terms have been omitted from Sq and our analytic compu- 
tations: the Hopf term which endows the Skyrmions with fractional statistics, and the long- 
range Coulomb interaction between the Skyrmions. We believe this is permissible because 
of the charge gap. Furtheiil, as the layers are antiferromagnetically correlated, Skyrmions 
in one layer will be correlated with anti-Skyrmions in the other, and this neutralizes the 
leading contribution of both terms. This latter argument should continue to hold even if 
the charge gap were to vanish at a quantum critical point (the charge gap remains non-zero 
at the quantum critical points in both our present calculations). Note also that no new 
term is necessary to induce charge transfer between the layers: a hedgehog/anti-hedgehog 
pair in the two layers corresponds to an event transferring Skyrmion number between them. 
Such spacetime singularities are absent in the strict continuum limit but appear when a 
short- distance regularization is introduced. 

For completeness, we note that the purely n field formulation becomes incomplete for 
m > 1 and larger g, as the spin zero Laughlin quasiparticles can become the lowest energy 
charged excitations. These should, in principle, be accounted for by a separate complex 
scalar field. However, these can also be neglected for the same reasons presented above for 
non-zero spin charged excitations. 
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We now manipulate the effective action into a form more suitable for our subsequent 
analysis. We solve the constraints 2 = 1 by representing 



(3.5) 



where n and L are vectors satisfying 



n 



L-n = (}. 



(3.6) 



Note that this representation is so far exact. Next, we insert ( p.5|) into Si. Because the 
layers are antiferromagnetically correlated we expect that L will not be too large, and it is 
therefore permissible to expand the resulting action to quadratic order in L. This is clearly 
an approximation: in Appendix |^ we examine a model solvable Hamiltonian by the same 
method in order to assess the damage done by neglecting terms higher order in L — we find 
that this procedure obtains the low energy spectrum correctly but introduces some spurious 
states at higher energies. To quadratic order in L, Si takes the form 



dr 



(3.7) 



where 5 is a unit vector in the direction of the magnetic field. Now we integrate out L 
while maintaining the constraint L ■ n = by adding an additional term to the energy 
C{L ■ nY and then taking the limit C —>■ 00. This yields the following effective action for 
the antiferromagnetic order parameter n 

\ 2" 



1S2 = — d X 
2t 



l/T 



dr 



1 / dfi 

(ySY + ^ ( ^ - i^zZ X n 



(3.8) 



where 



J 



1/2 



Mi 



(3.9) 



This is precisely the action of the 2+1 dimensional quantum 0(3) non-linear sigma model 
in a field B coupling to the conserved global 0(3) chargc^Hnl It is expected to apply to 
double-layer quantum Hall systems with z/ = 2/m at length scales larger than ~ do- 
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The remainder of this section consists of a detailed analysis of the properties of 1S2. The 
techniques and some results have already been presented earlier in Refs |2^J2^J2^: we shall 
present here a unified treatment with a special emphasis on dynamical properties at non-zero 
temperature. We begin in Section |111 A| by developing a simple mean-field phase diagram of 

^2. 



A. Mean field theory 

This section will summarize the results of the application of the mean field theory of 



Ref |27| to the action ^2. Formulation of the mean field theory requires some short distance 
regularization, and we choose to place the continuum theory on a square lattice in the spatial 
directions, with a lattice spacing a ~ a continuum formulation is maintained along the 
time direction. The resulting action is equivalent to the following lattice quantum rotor 
Hamiltonian 

H = ^ (^^L, - A J -t^-K J2 fi,- (3.10) 
where the couphng constants in 7i are 



<ij> 



f = -, K = ^ (3.11) 

t 

The Hamiltonian is expressed in terms of operators ni which represent the orientation of 
the rotors, and Lj which are the rotor angular momenta. The operators on different sites 
commute, while those on a single site obey the commutation relations (dropping the site 
indices) 

Lai Lf3 = iea/S'yLy 



[n„n^]=0 (3.12) 

We will describe the properties of Ti by choosing the best among the mean field Hamil- 
tonians given by0 
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Umf = E - ■ - ■ k)j (3.13) 

Here Z (= 4) is the lattice co-ordination number, and is a variational parameter to be 
chosen so that the expectation value of Ti in the ground state oH-Lmf is as low as possible; 
by the usual argument, this is expected to happen when N = (n). 



As in Ref ^ we numerically diagonalized Hmf by truncating its spectrum at some large 



angular momentum, and then optimized the value of N. The resulting phase diagram is 



shown in Fig 13. We discuss the properties of the various phases in turn: 



1. Spin Singlet (SS or SYM) 

Each rotor is in its non-degenerate i = state, = 0, and there is a gap to all excitations. 
The ground state is a spin singlet, and is therefore unaffected by variations in the value of 
A,. 

2. Quantized Ferromagnets ( QFi ) 

Again, A^ = 0, each rotor now has azimuthal angular momentum m = i and this value 
remains pinned as various parameters are varied. Each rotor is also in precisely the state 
with L = i{i +1), although this latter feature is a special property of mean field theory 
which will not survive fluctuations. Of these phases, only the i = 1 case is actually allowed 



for the double layer action Si, and it is clearly the FM phase of Fig |I0[ The other phases are 
an artifact of the approximations made in mapping Si to S2'. this should be clear from the 
discussion in Appendix ^ where we show that expanding in powers of L introduce spurious 
higher angular momenta states. 

3. Canted (C) and Neel (N) States 

These states have iV 7^ and varying continuously as the parameters are varied; we have 
Nx 0, Ny and A^^ = 0. From (|3^ ), this implies that the two layers have opposite spin 
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polarizations in the x — y plane. The two layers also have an identical ferromagnetic polar- 
ization, given by (L) which is oriented along the z direction. This ferromagnetic moment 
varies continuously as parameters are varied, and vanishes when ^ 0. So for general 
Az ^ this state is canted (C), while for = it is a pure Neel (N) antiferromagnet. 
The C phase has a single linearly dispersing spin wave mode in the x — y plane, while the 
N phase has two spin waves0. 

In the remainder of this section, we will present a detailed theory of the universal prop- 
erties of the system in the vicinity of the multicritical point M. This is the same quantum 



multicritical point (M) which exists in the universal Hartree-Fock phase diagram of Fig. |10 
where the N phase (along the A^ = line), the C phase and the SYM (SS) phase come 
together at Agas = 2?7_. We point out in this context that the other distinct multicritical 
point of the Hartree-Fock theory where the canted antiferromagnetic phase, the ferromag- 
netic phase, and the symmetric phase coexist (the point on the abscissa defined hj d = 
and A^ = Ag^s in Fig. ^) is not accessible within the effective field theory due to the long 
wavelength restriction d > Iq- (We mention that in our notations A^ in the Hartree-Fock 
theory corresponds to just B in the field theory due to our choice of units.) 

Note that the C, N, and SS phases meet at M, and so we will also discuss the universal 
second-order transitions between them. We will not discuss the nature of the second-order 
transitions between the QFi and C phases: very closely related transitions, in the same 



universality class, have been discussed in some detail in Ref 27 



B. Zero temperature critical properties 



A first study of the properties in the vicinity of the point M has appeared in Ref ^ 
using a large N expansion in a non-linear sigma model with component fields. The issues 
of interest here are more conveniently obtained using a recently developed expansion^ in 
spatial dimensionality, d, in powers of e = 3 — d. The latter approach is expressed in terms 
of a soft-spin field theory, and we therefore begin with a soft-spin version of the non-linear 
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sigma model 1S2: 

»1/T 



(3.14) 



Here (j) = {(px, 4'y, 4'z) ^ n is the soft-spin field which measures the staggered moment of the 
two layers. We have taken the magnetic field to point in the z direction. We will also be 
interested in the uniform ferromagnetic moment density of the system, M, and this is given 

by 

M = Moinu + n2.) = (3.15) 

where is the free energy density associated with the action S^p. We have introduced two 
new coupling constants, r and uq in Sip; these are related to the coupling t of 1S2, and its 
short-cutoff ~ io- We will not specify the precise values of these parameters here, as they 
merely appear at intermediate stages of our computation, and not in our final results. 

Let us first discuss the mean field properties of S^, obtained simply by minimizing the 
action while ignoring all spatial and time dependence of 0. For r — > 0, the ground state 
has (0) = 0, and is therefore in the quantum paramagnetic SS phase. For r — A^ < 0, the 
ground state has (0) 7^ and in the x — y plane. This is the C phase and the fields have the 
expectation values 

'^M^Y\o,o] M^5^i(M^. (3.16) 

Uo J J Uq 

or any rotation of in the x — y plane. Notice that M vanishes for A^ = 0, and therefore 
the line r < 0, A^ = is the N phase. The resulting mean field phase diagram is shown in 
Fig n. Notice that the vicinities of the points M are very similar in Figs [T^ and n. The 
quantum critical point M is at A^ = 0, r = 0, and it is clear from the Lorentz-invariant 
structure of at A^ = that this point has dynamic exponent z = 1. Rotations of the 
order parameter in the x — y plane have associated with them a stiffness ps(A^) given by 

Ps{A,) = (3.17) 
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This is the same stiffness which was computed in Section \Ll C\ and Eqn p.26| in the Hartree- 
Fock theory. 

We now include the effects of fluctuations at one loop. We will quote results for the 
dynamic longitudinal and transverse susceptibilities of the field which are measured in 
light scattering. Recall that in terms of the spin polarizations of the two layers ni , ^2 , we 
have ~ ni — ^2. We define (with T as the temperature ks = 1 in our units in this section) 

Xwiiu;) = J d^x fj^ dre-^'^^ (0,(a;, r)0,(O, 0)) (3.18) 

and 

Xdi^) = \j fj^ dre-'^^ ((0,(x, r) + i<j)y{x, r))(0,(O, 0) - z0,(O, 0))) (3.19) 



We can use the methods of Ref ^ to compute the one loop values of these susceptibilities 
in the SS phase (this is the phase with no broken symmetry) in the vicinity of the point M; 
we obtain 

Here the quantity A is an observable defined by 

A = Spin gap of the SS phase at T = for r > and A^ = 0. (3-21) 

The value of A should either be measured experimentally, or computed by a detailed mi- 
croscopic calculation like the Hartree-Fock theory discussed earlier in the paper. We will 
express all our results for r > completely in terms of universal functions of parameters A, 
T and B (so that the microscopic couplings r and uq do not appear anywhere in our results.) 
Clearly, in the mean field theory A = ^/r\ at one loop order, we have A ~ , where the 
exponent v = 1/2 + 5e/44. 

We need a separate experimental observable to measure the deviation of the system from 
the point M at A^ = for r < 0. A convenient choice, also used in Refs is the spin 

stiffness. We therefore define 
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Ps{0) = Renormalized spin stiffness of the N phase at T = for r < and = 0. (3.22) 

All our results for r < will be expressed in terms of Ps(0). Again Ps(0) ~ \r\'^, and the 
actual value of Ps(0) should be measured experimentally or computed in Hartree-Fock or 
microscopic numerical studies of the double-layer Hamiltonian. 

Before closing this subsection, we draw attention to the fact that there are two phase 
boundaries that terminate at the point M: the SS to C transition and the N to C transition. 
In the vicinity of these transitions the response functions computed near the critical point 
M should turn into reduced scaling functionsiiS characteristic of the respective phase tran- 
sitions. In the following subsections, we discuss simplified versions of the action iS,^ which 
can be used to compute these reduced scaling functions. 



1. SS-C transition, |A — A^l <^ A, r > 



In this region we can neglect (pz fluctuations and focus only on the 0^. + i(f)y which is 
undergoing Bose condensation. Further, it can also be shown that the second-order time 
derivative in can be dropped. Making these approximations, and defining 



(3.23) 



we see that reduces to 
S^, = j d^x J dr 



Uq 



9r 2AJ " ' ^ ' 24A2 



^1 



(3.24) 



This action has been previously studied in some detailE^a; j-^as a 2; = 2 quantum critical 
point at A = A2, and we will use the existing results later. Thus the SS-C transition is a 
line of z = 2 critical points terminating in 2; = 1 critical end-point M. 



2. N-C transition, B <^ ps{0), r < 

Both the N and C phases are ordered, and it is sufficient to simply focus on static, ther- 
mal, orientational fluctuations of the order parameter. We therefore quench the magnitude 
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fluctuations of (f) and return to tlie fixed length vector n. The effective action for static n 
fluctuations can be deduced from Sfj, to be 

Sn = ^l d^x fp.(0)(V,n)2 + 7^2] . (3.25) 



2T 

As noted earlier, Ps{0) is the spin stiffness of the Neel state, fully renormalized by quantum 
fluctuations. The anisotropy 7 = 6A^(A^ — r)/uo to lowest order in uq, and we expect 
7 ~ more generally. The action Sn has been studied in Ref and we will use their 
results in the following subsection. 

C. Non-zero temperature response functions 

A number of new phenomena occur at non-zero temperatures, and these are addressed 
in a relatively straightforward manner using the present continuum effective field theory 
approach. 

• There is a broken x — y symmetry in the C phase, and therefore a non-zero temperature 
(Tc) at which this order disappears in a Kosterlitz-Thouless transition. An estimates 
Tc was given earlier (sections p.1 C| and Fig. [Tl|) in the Hartree-Fock theory which is 
valid when the system is well away from one of the T = phase boundaries of the C 
phase in Figs O, and IT3. We expect to vanish continuously as the system in the 



C phase approaches the T = boundaries to the N or the SS phase: there is nonzero 
temperature phase transition above the N or the SS ground state. We discuss below 
the behavior of Tc near the C-N and C-SS T = phase boundaries. Near the point M, 
Tc is determined completely and universally by the two energy scales which measure 
the deviation of the ground state from M. So for r > we expect 

Tc = A,v^> ( (3.26) 

where \['> is a fully universal function; because the SS-C phase boundary occurs pre- 
cisely at A = A2, we have 
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^>(m > 1) = 0. (3.27) 

Similarly for r < we have 

T. = A.vI/< I^pM^ (3.28) 

where \E'< is also a universal function. Clearly the two functions should agree at r = 0, 
and therefore we have ^E'>(0) = \E'<(0); actually it is possible to say much more — for 
> we expect that is a smooth and analytic as a function of r through r = 0, 
and so using the dependencies of A and Ps(0) on r, it is possible to express \E'>,< 
as analytic continuations of each other. We will give explicit expressions for ^>,< to 
leading order in the expansion in e = 3 — d below. 

• The one-loop T = results for the SS phase ( |3.20| ) predict infinitely sharp absorption 
peaks in x\\ at u = A, and in at u; = A ± A^. As the SS phase has a spin gap, 
we expect these infinitely sharp peaks to survive at higher orders in the perturbation 
theory at T = 0. For T > two qualitatively new features will arise. First, thermal 
damping will lead to a broadening of the peaks. Second, the peak positions will 
themselves become temperature dependent. We will describe these processes below in 
the vicinity of the point M, where both the broadening and the T-dependent shifts 
are quite significant. Deep inside the SS phase, well away from the M point, these 
T-dependencies are exponentially activated, and therefore much weaker. 

We will restrict our results for the most part to the paramagnetic phase, although results 
in the magnetically ordered phases can be obtained by very similar methods. This means 
that we are working at T > above the SS phase, and at T > Tc above the C phase, all 
within the vicinity of the point M. The results are obtained using methods discussed in 



some detail in Ref ^ the only change is that the Zeeman splitting A^ has to be included in 
the propagators for the (f)x,y fields, and this modifies the values of the Matsubara frequency 
summations in the loop diagrams by replacing an energy £ by £ ± A^. The reader may also 
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consult Appendix D of Ref |3^ where a simpler derivation of just the one loop results of 
Ref ^ is given. 

The non-zero T generalization of (|3.20| ) takes the form 



2 



Here m\\^± and Fy^x depend implicitly upon the energy scales T, A^, and A (ps(0)) for r > 
(r < 0) in a manner we shall describe below to lowest order in e. Clearly, the "masses" ""^n.^ 
represent the peak absorption frequency, while Fy ^ are the absorptive pieces which lead to 
a T-dependent broadening of the line. 
First we describe the behavior of rn±^\\. 

For r > 0, the masses are universal functions A, T and A^. They can be written as 

2nT 
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(3.30) 



The function G{y, h) represents the value of the one-loop momentum integral; it was com- 
puted in Refs for the zero magnetic field case h = 0. The generalization to non-zero 
h is 



G{y,h) 



-2 / dq 
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(3.31) 



q^ + y-h^ ) - 2^q^ + l/e_ 

This integral has to be evaluated numerically in general, but we have the limiting value 
G(0,0) = 27r^/3. Stability of the paramagnetic state requires that m± > A^; this require- 
ment leads to an expression for T^, which is determined by solving m± = A^. Analysis of 
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this equation in powers of e shows that Tc ~ This imphes that A/T, A^/T ~ y/e, and 

so to leading order we can just use the value of G{0, 0) in ( p.30|) to obtain 



_ 33(Ag - A^) .3 

for Az > A. For A^ < A the system is in the SS phase, and therefore = 0. Notice that 
( p.32| ) agrees with the scaling form ( p.26| ). This result is expected to be the leading order 
result in powers of e, except in the region | A^ — A| ^ A where the e expansion fails and the 



reduced action Sq, of Section |1I1 B 1| has to be used. Using results of Ref ^ for the latter 



action we have the exact asymptotic form 
(A,- A)ln[A,/(A,- A)] 

Closely related results can be obtained for r < 0. In this case, the masses are universal 
functions of Ps(0), A^ and T. However, considerable ambiguity arises in the e expansion 
for the result because Ps(0) does not simply have the dimensions of energy for all d. The 
appropriate scaling variable^ is {ps{0)Y^^'^~^\ and it is necessary to keep the full l/{d — 1) 
power, rather than expand it in powers of e in order to make the engineering dimensions of 
the results come our correct. This then leads to ambiguities as to precisely which numerical 
factors should be raised to the power l/{d — l) and which to 1/2 + e/4. A convenient choice, 
which leads to the most compact expressions is to define 

^= r4^^ (3.34) 
\{n + 8) Sd+i J 

where we have written the general expression for the n-component order parameter: in 
the present case n = 3. The factor Sd+i is a phase-space factor and is given by Sd = 
2/[r(rf/2)(47r)"'/2] (this factor was inadvertently omitted in Ref ^). Notice that ps has the 
dimensions of energy in = 2 (which is of interest here). The value of p^, however must be 
regarded as subject to large systematic corrections, in view of the ambiguities noted above. 
Using the methods and results of Ref U^ for r < 0, and expressing them in terms of p^, we 



find that the results ( p.29|) still hold, but (|3.30|) are replaced by 
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(3.35) 



Notice that G{y, h) is now needed for negative values of y. Despite appearances, the expres- 
sion ( ^.31[ ) actually also holds for y < — one simply uses the identity cosh(zx) = cos(a;) 
when the square root becomes purely imaginary. Indeed, it is not difficult to show that the 
expression in (|3.31|) is actually analytic for all real — oo < y < oo provided h > 0. We can 
use the same stability condition used for r > to now obtain the leading order e-expansion 
result for Tc. 

33{Al + pi/2) 



(3.36) 



which is of the scaling form (|3.28|) . The e expansion fails when ^ Ps(0) where the system 
approaches the C-N phase boundary; here, we use the effective action Sn of Section [lllB 2 
and results for it in Ref ^ to obtain 

2vrp,(0) 



for ln(p,(0)/A,) > 1. 



(3.37) 



ln(p.(0)/A,) 

Finally, we obtain the damping coefficients ^±,\\- This requires evaluation of two-loop 
diagrams and the results are extremely lengthy. We will be satisfied here by simply quoting 
the results valid for A^/T ^ 1, (A or ps{0))/T ^ 1 which were obtained in Ref 



r,M = riiM 



121 



UJ 



^ + 7r2T2 + 6T2Li2(e--/2T) 



(3.38) 



where Li2(a;) is the dilogarithm function 

U,ix) = - r^Hl-y) 
Jo y 



(3.39) 



D. Connection to the Hartree-Fock theory 



The effective field theory for the double-layer QH system at a filling factor oi u = 2/m 
( with m an odd integer) that we develop above is entirely built on the effective action 

40 



5*2, defined by Eq. (|3.8| ). In particular, we make use of tlie fact that this effective action 
for our problem is identical to the action of the 2 + 1 dimensional 0(3) quantum non- 
linear a modelS^ with the additional feature of an external magnetic field coupled to the 
conserved global 0(3) charge. Once this precise mapping of our effective action to that 
of the 2 + 1 dimensional 0(3) quantum nonlinear a model becomes explicit, the rest of 
the results derived in sections |l|A-C follow naturally. The question now arises about the 
correspondence between our effective field theory results in this section and the microscopic 
Hartree-Fock results (for v = 2) described in sections [II AHII C^ . 



It is to be noted that both the microscopic Hartree-Fock theory (sections [II A| - |II C]) and 



the effective nonlinear a model field theory predict the same number of zero temperature 
quantum phases, namely the fully spin polarized ferromagnetic, the canted antiferromag- 
netic, the Neel, and the symmetric spin singlet phase, for the double-layer QH system at 
V = 2. (The effective field theory, in addition, enables us to predict that the double-layer 
system at all fillings v = 2z/i, where z/i = l/m with m odd, has these four phases with the 
spin singlet phase in the general case being a non-trivial correlated SS phase rather than 
just the pseudospin-symmetric spin-antisymmetric SYM phase of the u = 2 Hartree-Fock 
theory.) It should also be noted that both the Hartree-Fock theory and the effective field 
theory predict the existence of a finite temperature Kosterlitz-Thouless phase transition in 
the canted antiferromagnetic phase with the planar antiferromagnetic ordering disappearing 
above the Kosterlitz-Thouless transition temperature. The underlying physics in both the 
theories is that the system is essentially an X-Y antiferromagnet in the layer (i.e. in the 
plane perpendicular to the magnetic field direction) in the new canted phase. 

On a more quantitative level it is easy to show that both theories predict the same 
topology of the zero temperature quantum phase diagram. This is demonstrated in Fig. 
p!5[ where we have redrawn the Hartree-Fock phase diagram (Fig. ^a) of Fig. [1^ inverting 
abscissa (from Agas/U- to U-/Asas) and have somewhat reconfigured the effective field 
theory phase diagram (Fig. ^^) from Fig. by keeping only the QFi phase and by 
modifying the relative size of the various phases (which are arbitrary within the effective 
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field theory). Using the definitions t = {J/ p^M^Y^"^ from Eq. ( p.9|) to define the abscissa, the 
effective field theory phase diagram in t — space (Fig. p!5|b) can be seen to be identical 
topologically to the quantitatively calculated Hartree-Fock phase diagram (for u = 2) in 
the Asas/U- — Az/Asas space (Fig. 0). Note that, in addition to the identical topology 
involving four distinct quantum phases as shown in Fig. ^ of the two phase diagrams 
with the effective coupling parameter t of the field theory (the abscissa in Fig. p!5|b) being 
proportional to the parameter A^as/U- (the abscissa in Fig. p!5|a) of the Hartree-Fock theory 
(which is expected, because t ~ Agas/U with J being the interlayer superexchange couphng) 



and the ordinate (~ A^) being the same in both Figs. |T5]a andp!5|b, the multicritical point M 
on the zero magnetic field line shows up in both phase diagrams. At the (zero temperature) 
quantum multicritical point M, the canted, the spin singlet, and the Neel phase coexist. 
(The other distinct multicritical point of the Hartree-Fock theory, which is apparent on the 
abscissa of Fig. |l^ where A^as = A^ and d = 0, where the ferromagnetic, the canted and 
the symmetric phase coexist is not accessible within the effective field theory because of its 
long wavelength approximation, and can not be seen in Fig. pTS^ as it is pushed to the point 
U-/ Asas = 0, A^/Asas = 1 whcre the two Hartree-Fock phase boundaries of Fig. p!5| a come 
together.) It is, therefore, obvious that, except for very small values of d (where the effective 
field theory which applies only when d > lo), the quantum phase diagrams predicted by the 
two theories are topologically identical. 

Finally, we can actually estimate the u = 2 double-layer Kosterlitz-Thouless transition 
temperature, of section pT[|C, in the effective field theory by using the microscopic param- 
eters obtained within the Hartree-Fock theory. This calculationi, where one incorporate the 
calculated Hartree-Fock parameters for A and e = 1 in Eq. O), leadsil to an estimated 
effective field theory Tc ~ 3K which is somewhat larger than the critical temperature Tc (Eq. 
( p.27| )) estimated within the long wavelength mean field Hartree-Fock treatment of section 
P C| . In general, we believe the e-expansion leads substantial overestimates of transition 
temperatures because it does not properly account for the low- dimensional vortex effects 
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responsible for the transition. 



IV. COMPARISON WITH EXPERIMENTS 

In this section, we discuss some recent double-layer v = 2 inelastic light scattering ex- 
periments whose findings are consistent with our theoretical results. A detailed quantitative 
comparison between our theory and the experiment requires an accurate knowledge of the 
temperature dependence of the related experimentally relevant response functions as the sys- 
tem undergoes a finite temperature phase transition at T^. Such a quantitative description 
is at present lacking, and therefore we restrict ourselves mostly to a qualitative discussion. 

In a recent inelastic light scattering experiment,! the long wavelength uq mode of the 
intersubband SDW triplet (see Fig. |T6| for schematic details of the various possible SDW 
modes in the system), which corresponds to transition |Ocr) is measured for v = 2 

double-layer quantum Hall systems. The double-layer samples used in the experiment are by 
design in the canted antiferromagnetic phase according to our zero temperature Hartree-Fock 
phase diagram, i.e., the ground state of the experimental system is the canted antiferromag- 
netic quantum Hall state (see Figs. |^, |10|, and |15| for the location of the experimental sample 



in our theoretical diagram). The experiment^ shows two important and striking features: 
One is that there is a threshold temperature (~ Q.bK) below which the ujq mode becomes 
unobservable as it seems to lose all spectral weight, the other feature is that the excitation 
energy ujq approaches the Zeeman energy when the threshold temperature is approached 
from the above, i.e. ujq ~ A^. We argue below that these experimental observations are 
completely consistent with our predicted Kosterlitz-Thouless transition in the canted anti- 
ferromagnetic phase being the observed experimental transition at Tc. 

First, we notice that the ujq mode, which involves a no-spin-fiip transition with 5Sz = 0, 
has a maximum spectral weight in the symmetric phase, where there are as many spin-up 
(down) empty states as there are spin-up (down) electrons. The spectral weight of the ljq 
mode is identically zero in the ferromagnetic phase, where all spin-up states are occupied 
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and all empty states are spin-down, and hence the ujq mode (which does not involve any spin 
flip) is forbidden. The spectral weight of the ujq mode should be nonzero but small in the 
antiferromagnetic phase. This is because the canted antiferromagnetic phase lies between 
the symmetric phase and the ferromagnetic phase in the phase diagram and its spin-flip 
dynamics should thus be intermediate. More over, the canted antiferromagnetic phase is 
not an eigenstate for either spin or pseudospin, so the small spectral weight of the ujq mode 
is shared by many allowed transitions, spreading the mode intensity over these transitions 
and thus making the spectral weight of each transition even smaller. It is thus plausible 
to regard the observed disappearance of the uq mode at the threshold temperature as the 
transition to the canted antiferromagnetic phase at lower temperatures (where the spectral 
intensity for the Uq mode becomes very small). Above the transition temperature the system 
is essentially a disordered planar X-Y magnet, and thus behaves like a paramagnet whose 
SDW properties should be very similar to the paramagnetic spin-singlet symmetric phase. 

Next, we notice that, in the symmetric phase, the excitation energies of the intersubband 
SDW triplet have the following simple relationship 

cj± = cuo±A,. (4.1) 

This expression can be derived explicitly, using either the diagrammatic time- dependent 
Hartree-Fock approximation or the single-mode approximation. It is a direct consequence 
of the fact that Coulomb interaction is spin independent. The above relationship bears a 
clear physical meaning: — ^ means that uj_ 0, i.e. mode softening (see Fig. |16|). 



Thus, the experimental observation that Uq approaches the Zeeman energy as the threshold 
temperature is reached from above suggests that there is mode softening {u^ = 0) at the 
phase boundary, as predicted by the computations of the T-dependent peak positions in 
Section PK^l 

Finally, we note that the critical temperature (the threshold temperature) in the 
experiment! is ~ 0.52K, which is reasonably close to our calculated Kosterlitz-Thouless 
critical temperature Tc ~ 1.8K in the Hartree-Fock theory (Eq. ( p.27|) ) using the ac- 
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tual experimental sample parameters. □ This discrepancy between the experiment and the 
Hartree-Fock theory is small when compared with the energy scale of Coulomb interaction, 
which is about 707^ in this particular sample. In addition, quantum fluctuations neglected 
in the Hartree-Fock theory should lower the calculated critical temperature and reduce this 
discrepancy. 

From the above discussions, we conclude that our theoretical predictions are consistent 
with the recent light scattering experimental results. The most dramatic aspect of the 
experimental observations which give us confldence in believing that the experiment is really 
seeing the canted antiferromagnetic phase are (i) the unambiguous observation of a mode 
softening {i.e. implying — > 0); (ii) the observed temperature dependence 

indicating a flnite temperature phase transition; (iii) the location of the experimental sample 
in our calculated phase diagram and (iv) the uq — > collapse being observed precisely at 
z/ = 2. 

While the recent inelastic light scattering experimenti provides, in our opinion, rather 
compelling evidence in favor of there being a flnite temperature (Kosterlitz-Thouless) tran- 
sition in the u = 2 double-layer system with the low temperature phase being the canted 
antiferromagnetic phase (by virtue of the vanishing of the tu. mode at the phase boundary, as 
discussed in section [II B| of this article), a complete veriflcation of our theory awaits further 
more conclusive and direct experimental measurements, especially heat capacity measure- 
ments which should shown (Fig. |12D power law temperature dependence in the canted phase 
due to the existence of the Goldstone mode and exponential temperature dependence in the 
two normal phases due to the existence of gaps in the excitation spectra, would be par- 
ticularly well-suited in verifying our phase diagram. The direct observation of a gapless 
Goldstone mode (Fig. |^) in the inelastic light scattering measurement in the (low tempera- 
ture) canted phase would also be rather deflnitive in establishing the existence of the canted 
phase. In this context we mention that the SDW softening indicating a phase transition to 
the canted phase is a long wavelength instability, and therefore optical spectroscopy^ may 
also be useful in studying our proposed u = 2 double-layer phase diagram. Both of these pro- 
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posed direct experiments are fraught with considerable (experimental) difficulties, however. 
Electronic heat capacity measurements in quantum well structures are notoriously difficult 
by virtue of the extremely small magnitude of the (2D) electronic heat capacity compared 
with the background (lattice) contribution. As for the direct experimental observation of 
the Goldstone mode, the experimental inelastic light scattering spectroscopy is severely re- 
stricted by the selection rules inherent in the resonant light scattering spectroscopy, and at 
this stage it is unclear whether the problems associated with the selection rules would allow 
to directly observe the Goldstone mode. 

One striking difference between the physics of = 2 double-layer system and the cor- 
responding u = 1 situation is the existence of a charge gap in the u = 2 case for all values 
of d and Agas^ the system is always incompressible (in all its quantum phases including the 
canted phase). Thus the quantized Hall effect exists throughout our phase diagram unlike 
in the corresponding u = 1 situation.i00'Ei The existence/nonexistence of the QH effect, 
which has been useful in mapping out the u = 1 double-layer phase diagram^ would not 
work in our problem in a direct sense. We do, however, speculate that the activation en- 
ergy {i.e. the effective value of the incompressible charge gap) for the u = 2vi double- layer 
QH effect may very well show observable structure at our calculated phase boundaries even 
though all the phases (ferromagnetic, canted, symmetric) would exhibit v = 2vi QH effect. 
We suggest systematic experimental investigations of u = 2vi double-layer { vi = 1/m with 
m = 1,3,5, ...) QH activation energies by tuning Agas, A^, and d to look for signatures of 
our proposed zero and finite temperature phase transitions. 

In this context we point out that there is already some experimental evidenceS'll that 
the naive A^ = Agas level crossing in z/ = 2 double-layer QH systems does not exist (as 
our theory proposes and clearly demonstrates in our calculated phase diagrams). The ex- 
perimental observationii'ii has been that the naive v = 2 level crossing phenomenon (at 
Az = Asas) between ferromagnetic and symmetric phases, which would exist in the absence 
of our intervening canted phase, if it happens at all in double-layer systems, must happen 
at magnetic fields much lower than that satisfying A^ = Agas condition. This is, of course, 
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exactly what our phase diagram (see Fig. ^)) predicts — nothing interesting happens for 
finite d at Az = Agas or for that matter even for = SA^ at d = 21^ in Fig. ^a for 
example — the system remains in the fully spin polarized ferromagnetic phase and the naive 
expectation of a level crossing transition to the symmetric phase simply does not occur. 
In this sense, our phase diagram for the v = 2 double-layer system may have already been 
verified in 1990!i Further experiments along this line at u = 2z/i double-layer systems would 
be useful. 

V. SUMMARY 

In summary, we have studied both zero and finite temperature properties of the u = 2 
double-layer QH systems within the framework of Hartree-Fock approximation. We show 
that, in addition to the fully polarized state adiabatically connected to the well separated 
layer state, there are two other double-layer quantum Hall phases: the first is a spin singlet, 
and the second is characterized by a finite interlayer inplane canted antiferromagnetic spin 
ordering. The transition between the different quantum Hall phases is continuous, and is 
signaled by the softening of collective intersubband spin density excitations. Because of 
the broken U{1) symmetry in the canted antiferromagnetic phase, the system has a finite 
temperature Kosterlitz-Thouless transition (T^ ~ IK). Below the critical temperature, the 
canted antiferromagnetic phase supports a linear Goldstone mode. Above, the system is 
essentially a paramagnet similar to the symmetric phase. Our findings are consistent with 
recent light scattering spectroscopic experimental results. We present detailed results of 
our study, including the stability energetics of various phases, the antiferromagnetic order 
parameter in the canted phase, the phase diagram, the collective excitation dispersions, the 
specific heat, and the Kosterlitz-Thouless critical temperature, and suggest various experi- 
ments which could, in principle, probe the rich double-layer phase diagram predicted by our 
theory. 

In addition to the microscopic u = 2 Hartree-Fock theory, we have developed a rather 
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general long wavelength effective field theory for the v = 2z/i, where z/i = 1/m with m 
an odd integer, double-layer system. The essential inputs for this effective field theory are 
the existence of charge gaps in the two layers and an effective interlayer antiferromagnetic 
(superexchange) interaction. By mapping the effective action for this problem to that of 
an 0(3) quantum nonlinear a model, we have been able to show that the qualitative phase 
diagram calculated within the Hartree-Fock theory for u = 2 is actually generically valid 
(topologically) for any z/ = 2z/i (with z/i = 1,1/3,1/5,...) double-layer system with the 
symmetric phase of the Hartree-Fock calculation being replaced by a highly non-trivial cor- 
related spin singlet phase ( of which the u = 2 symmetric phase is a rather trivial example). 
Thus, there could be rather non-trivial canted (and perhaps even Neel, if one can apply suf- 
ficient external pressure to produce vanishing gyromagnetic ratio) antiferromagnets at, for 
example, z/ = 2/3 in a double-layer system, where each single fully spin polarized Laughlin 
state spontaneously develops in-plane antiferromagnetic spin ordering. Observation of the 
canted or the spin-singlet phase in a = 2/3 double-layer QH system would significantly 
enrich the many-body strong correlation physics associated with QH systems. 

We conclude by pointing out that, although we have confined ourselves in this article to 
the V = 2/m case, with m an odd integer, it is obvious that the physics we are considering 
here applies, in principle, to all double-layer QH systems with v = 2vi where a single layer at 
filling Vi forms a fully spin polarized incompressible QH state with a charge gap. Thus, the 
same physics as at z/ = 2 should apply, in principle, at z/ = 6 (but not at z/ = 4, 8, ... where 
the charge gap is the cyclotron gap not A2,Asas-) in double-layer system. In principle, 
however, our approximations which neglect all (orbital) Landau level couplings become 
progressively worse at higher Landau levels. In this respect, it is very gratifying that the 
experimental light scattering measurementsi find qualitatively similar (but quantitatively 
much suppressed) behavior at z/ = 6 as at z/ = 2, but the z/ = 4 situation is qualitatively 
different. 
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APPENDIX A: TWO SPIN PROBLEM 



Here we will assess the validity of the mapping from the action Si in to S2 in (^]8| 
by examining a simple toy model of two spins. We consider the Hamiltonian 



n = JSi-S2- A J ■ {Si + ^2) 

where S^i 2 are two quantum spins of spin S. The energy spectrum of B is clearly 
J 



(Al) 



Ee = + 1) - A,m + Eq 



0,1,... 25 ; m = + 



- 1, 



(A2) 



where Eq is an overall constant we shall not be interested in. 

Let us attempt to obtain this result using the coherent state path integral. First, we 
transcribe Ti. into the effective action 



S 



dr 



iSA{ni) ■ ^ + iSA{n2) ■ ^ + JS^ni ■ 



(A3) 



dr dr 

where 2 = 1- Notice that this is the analog of the action Si in (|3.4|) with only the spatial 
gradient spin stiffness terms now being absent. Now insert the parameterization ( p.5|) into 
(|A3|), and expand to quadratic order in L. The neglect of terms higher order in L is the 
only approximation being made here. This gives us the analog of (|3.7|) 



S 



dr 



-> / dn 
2iSL ■ {^^'q — ^ ^^2-^ 



(A4) 



Now we integrate out L as described above ( p. 81) to obtain 



(A5) 
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where recall that the functional integral is over the unit vector field n{T) satisfying fp = 1 
for all r. This last form of S is the effective action for a quantum rotor in a field A^z. This 
action is equivalent to the Hamiltonian 

Hr = -AJ-E (A6) 

where L is the rotor angular momentum operator. The eigenvalues of T-Cr are easily seen to 
be identical to those of Ti in (|A2|) with one simple difference. The allowed values of i now 
extend over all non-negative integers. Thus the only effect of dropping terms higher order 
in L in the functional analysis is that the upper bound i < 2S has disappeared. This only 
introduces additional states at relatively high energies and is therefore not expected to be 
of importance in our study of the low energy properties of S2. 
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FIGURES 

FIG. 1. The energy per magnetic flux in the symmetric (SYM) state, the spin polarized 
ferromagnetic (FM) state, and the canted antiferromagnetic (C) state for a z/ = 2 double-layer 
system with Asas = 0.07e^/e/o, = O.Ole^ /elo, and the well-thickness dw = O.Slo- 

FIG. 2. Schematic display of electron spin orientations in the canted antiferromagnetic quan- 
tum Hall phase. 

FIG. 3. The canted antiferromagnetic order parameter versus layer separation for the indicated 
tunneling and Zeeman energies. The well-thickness d = O.Slo- 

FIG. 4. The zero temperature phase diagrams a,tv = 2 within the Hartree-Fock approximation 
at two different values of the Zeeman energy: (a) A^ = O.OleVe^o and (b) A^ = O.OleVe^o- The 
quantum well thickness is dyj = O.S/q for both the figures. Three phases are present: the symmetric 
phase (SYM), the spin polarized ferromagnetic phase (FM), and the canted antiferromagnetic phase 
(C). The '-|-' in (a) denotes the experimental sample parameters of Ref. |8[ The vertical dotted line 
in each figure indicated the A^ = Asas condition, which is the naive phase boundary between the 
FM (A^ 

^ Agas) a,nd the SYM (A^ < Agas) phase with an expected level crossing at A^ — Agas^ 

FIG. 5. Feynman diagram for the intersubband spin density response function in the 
time-dependent Hartree-Fock approximation, where solid lines are the self-consistent Hartree-Fock 
electron Greens functions and wiggled lines are Coulomb interaction potentials. 

FIG. 6. The inter-subband SDW dispersion a;±(g) in the spin polarized ferromagnetic (FM) 
phase at = 2 with tunneling energy Agas = 0.02e^/e/o, Zeeman energy A^ = O.Ole^/e/o, layer 
separation d = 1.15/o, and the well-thickness dyj = O.S/q. 

FIG. 7. The inter-subband SDW dispersion uj±{q) in the symmetric (SYM) phase at = 2 with 
layer separation d = 0.85/o, Zeeman energy A^ = 0.08e^/e/o, tunneling energy Asag = 0.35e^/e/o, 
and the well-thickness d^ = 0.8/o- 
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FIG. 8. The low energy intersubband SDW mode (^-{q = 0) and the canted antiferromagnetic 
order parameter (COP) versus tunneling energy with layer separation d = l.OZo, Zeeman energy 
Az = 0.08e^/e/o, and the well-thickness = O.Slo- 

FIG. 9. The inter-subband SDW dispersion w±((?) in the canted antiferromagnetic (C) phase 
at u = 2 with layer separation d = 1.15/o, tunneling energy Agas = 0.14e^/e/o, Zeeman energy 
Az = O.Ole^/e/o, and the well-thickness d^ = O.SIq- 

FIG. 10. Zero temperature phase diagram oiau = 2 double-layer quantum Hall system within 
the Hartree-Fock approximation. The phase diagram is expressed in terms of scaled dimensionless 
variables. The mark represents the experimental sample of Ref. |8[ The Neel phase (N) at 

= and Agas < 2U- is represented by the thick line. The M-point represents the quantum 
critical point at A^ = 0. 

FIG. 11. The calculated Kosterlitz-Thouless critical temperature Tc versus tunneling energy 
Agas at different interlayer separations: dotted line d = 1.4 /q, solid line d = 1.2 Iq, and dashed line 
d = 1.0 lo- Zeeman energy A^ = 0.04 e'^/elo- The layer-thickness d^ = 0.8 lo- 

FIG. 12. The heat capacity per magnetic flux of a u = 2 double-layer quantum Hall sys- 
tem as functions of temperature in the symmetry phase (SYM), in the spin-polarized ferromag- 
netic phase (FM), and in the canted antiferromagnetic phase (C). The inset shows C/T^, where 
T = kBT/{e'^ /elo), versus T in the C phase. 

FIG. 13. Mean field phase diagram of the quantum rotor Hamiltonian TC in ( p.lO| ). The phases 



are described in Section [III A| . Only the QF^ phase is expected to appear for the two-layer model 
under consideration here, and is referred to elsewhere as the FM: the QF2 phase is an artifact of 
the approximations made in deriving the rotor model. The SS phase was also called the SYM 
phase in the Hartree Fock computations. 
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FIG. 14. Mean field phase diagram of the soft-spin action 5^ in ( 3.14 ). The SS phase was also 



called the SYM phase in the Hartree Fock computations. Notice that it captures the vicinity of the 



point M in the rotor mean field phase diagram in Fig 13. The multicritical point M is described by 
a relativistic continuum field theory with dynamic exponent z = 1. The SS-C boundary is a line 
of second order transitions with dynamic exponent z = 2 and is described by action cSiji in ( p. 24 ). 



The position of this boundary is given exactly by A = A^, where A ~ r"^ is the A^ = spin gap of 
the SS phase {u is the correlation length exponent of M). The action Sqj holds for |A — A^l ^ A. 
The N state has T = spin stiffness Ps{0) ~ {—i^Y^ and for A^ ^ Ps{0), the action Sn in ( p. 25 ) 
describes low T fluctuations. 

FIG. 15. (a) The Zero temperature phase diagram of a double-layer quantum Hall system 



at u = 2 within the Hartree-Fock approximation. This is the same diagram as Fig. |10|. It is 
redrawn here with the abscissa inverted. The '+' mark represents the experimental sample of Ref. 
p. The Neel phase (N) at A^ = and Agas < 2U^ is represented by the thick line, (b) Zero 
temperature phase diagram of a double-layer quantum Hall system at v = 2ui derived from the 
effective Lagrangean 52 (Eq. ( |3.8D ). The inset shows the topologically identical Hartree-Fock phase 
diagram of Fig. |l^. The FPF, C, and SS phases in the main figure correspond to the FM, AF, and 
SYM phases in the inset, respectively. 

FIG. 16. The intersubband spin excitation transitions in a double-layer quantum Hall system 
at V = 2\n the (a) symmetric phase, (b) ferromagnetic phase, and (c) the canted antiferromagnetic 
phase. The spin conserved transition [loq mode) has large spectral weight in the symmetric phase 
and is prohibited in the ferromagnetic phase. 
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